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CHAPTER  1 


INTRODUCTION 


The  application  of  charged  particle  beam  drivers  to  inertial  confinement 

fusion  has  spurred  interest  in  the  study  of  energy  deposition  by  intense  beams 

in  extreme  states  of  matter.  Enhancement  of  energy  deposition  by  beams  has  been 

predicted  on  the  basis  of  various  modifications  of  single-particle  energy 

deposition  theory,  for  example:  collective  beam-target  interactions,  finite 

target  temperature  effects,  and  beam  stagnation  in  the  target.  A  review  of 

1  2  3 

these 'topics  has  been  given  by  Toepfer.  Observations  by  Imasaki  et  al.  * 
have  shown  that  enhanced  deposition  by  electrons  occurs  in  thin  targets  and  that 
there  is  a  strong  dependence  oY-the  effect  on  target  material.  They  have 
suggested  that  the  observations  are  explained  by  the  collective  interaction  of 

the  beam  with  the  target.  Recently  an  experiment  and  calculation  by  Young  et 

A  0  «  • 

al.  have  demonstrated  enhanced  energy  deposition  by  a  deuteron  beam  in  thin 

targets.  The  calculation  indicates  that  the  enhances  it  is  caused  by  the  free 
target  electrons  which  are  ionized  by  the  early  part  of  t’  aam  pulse.  Another 
example  of  increased  energy  deposition  is  the  observation  by  Kiselev  et  al.^ 
of  the  range  shortening  for  a  2  MeV,  1A  beam  in  air.  They  proposed  that  a 
collective  interaction  of  the  beam  with  the  plasma  that  it  generates  could 
explain  this  result. 
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McCorkle  and  Iafrate^  have  suggested  that  single-particle  energy 

deposition  theory  should  be  augmented  with  a  beam-density  term  when  intense 

beams  interact  with  a  plasma.  The  beam-density  effect  has  its  origin  in  the 

theory  of  stopping  power  of  molecular  ion  clusters,  given  by  Brandt,  Rathowski, 
7-9 

and  Ritchie,  which  is  in  turn  closely  connected  with  wakes  induced  in  the 

conduction  electrons  of  solids.^  McCorkle  and  Iafrate^  modified  the 

nonrelativistic  molecular  cluster  theory^  so  that  it  contained  a  characteristic 

relativistic  screening  length  and  used  the  result  to  estimate  the  beam-density 

contribution  to  energy  loss  by  relativistic  beams  in  a  plasma.  They  found  that 

the  effect  could  be  quite  large  relative  to  single-particle  energy  loss. 

13 

McCorkle  and  Cox  have  investigated  the  relationship  between  these  beam 
concepts  and  ohmic  loss  and  return  currents.^  ^ 

A  detailed  derivation  of  energy  loss  by  clusters  of  relativistic  particles 

19  20 

has  been  given  by  Rule  and  Cha,  which  includes  the  Fermi  density  effect,  and 

which  reduces  to  the  nonrelativistic  theory  in  the  appropriate  limit.  Recently 
21 

Rule  and  Crawford  have  shown  that  the  nonrelativistic  energy  loss  expression 
for  a  cluster  contains  long-range  dipole-like  terms  and  that  this  results  in  a 
beam-density  effect  which  depends  on  the  beam  shape  in  an  essential  way. 

The  present  paper  contains  a  detailed  theory  of  the  beam-density  effect  for 
finite  relativistic  beams,  starting  from  the  results  of  Reference  19.  An 
example  based  on  this  theory  is  given  for  a  50  MeV,  10  KA  pulse  of  electrons 
interacting  with  weakly  ionised  nitrogen.  It  is  shown  that  the  beam-density 
effect  results  in  a  large  enhancement  of  energy  loss  for  the  regime  of  beam  and 
medium  parameters  studied.  Finally  a  brief  discussion  of  the  relationship  of 
this  effect  to  other  related  phenomena  is  given. 
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CHAPTER  2 

CALCULATION  OF  THE  BEAM-DENSITY  EFFECT  FOR  A  FINITE  BEAM  PULSE 


THE  INTEGRATION  OVER  THE  BEAM  DISTRIBUTION  FUNCTION 


In  Reference  19  a  general  formula  for  the  coherent,  cooperative  energy  loss 


by  a  relativistic  pair  of  particles  was  derived.  The  case  of  "v^  parallel  to 


v^,  but  v.  >*  Vj  was  given,  where  v^  and  v^  are  the  velocities  of  the  particles. 
For  the  special  case  in  which  "v.  -  V.  **  V,  a  simpler  result  was  obtained  (see 
Equation  (3.33)  of  Reference  19),  namely 


(1-B2e (v)) 


(2.1) 


2  2 

The  quantities  Z.e  and  Z ,e  are  the  charges  of  the  particles;  u  ■  Avne  /m, 

1  J  P 

where  n  is  the  total  number  of  medium  electrons  per  unit  volume,  both  hound  and 
free,  and  m  is  the  electron's  mass;  v  is  the  frequency  in  units  of  c(v) 

is  the  dielectric  function  of  the  medium;  and  B  *  v/c.  The  coordinates  *  and  b 
are  the  magnitudes  of  the  components  of  the  vector  R^j  from  particle  i  to 
particle  j,  such  that  z  is  parallel  to  7  and  b  is  perpendicular  to *v .  The 
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quantity  "a"  in  (2.1)  is  the  minimum  impact  parameter  which  may  be  suitably 

chosen  to  take  into  account,  to  a  good  approximation,  the  close  collision 

portion  of  the  energy  loss.  The  modified  Bessel  function  Rq(Z)  appears  in 

(2.1).  For  the  special  case  of  particles  traveling  in  a  line  parallel  to  V,  an 

20 

analytic  expression  for  (2.1)  has  been  obtained. 

The  limiting  case  in  which  z-*0  and  b+0  in  (2.1)  gives 

W(0,0)  -  2  ZjZjS,  (2.2) 


where 


S  a  1  e2 

IT 


(?)2  Re  /  idv  v  (ifo  '  °2)  (1-s2‘)I/2)' 


(2.3) 


This  is  the  usual  expression  for  the  ionization  energy  loss  per  unit  length  of  a 

22 

relativistic  particle  of  unit  charge,  including  the  Fermi  density  effect. 


Using  (2.1)  and  (2.3),  one  can  write  the  expression  for  the  total 
ionization  energy  loss  of  a  pair  of  particles  separated  by  the  distance  as 


W  (?..)■  Z?S  +  Z*S  +  W(z,b) . 

-  ij  l  j 


(2.4) 


In  the  limit  R.j-*0  one  obtains,  using  (2.2), 


2  2 

lim  W(R£j)  -  ZiS  +  ZjS  2Z£ZjS  ■  (Z*  ♦  Zj)2S, 


(2.5) 


RirO 
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i.e.  the  energy  loss  for  s  particle  of  charge  ( Z ^  ♦  ZJe. 


With  reference  Co  (2.4),  Che  coCal  energy  loss  of  particle  i  can  be  written  as 


W.  =  Z^S  +  4w(z,b),  (2.6) 

— l  l  2 


and  similarly  for  particle  j. 


If,  instead  of  two  particles,  a  beam  of  particles  is  considers  t..en  the 

total  energy  loss  per  unit  length  by  the  ith  particle  in  the  beam  c  i  be 
written,  in  analogy  with  (2.6),  as 


-B 


Z?S 

l 


w_ 


(2.7) 


where 


U 

B 


dbb  p (z,b)W(z,b)  . 


(2.8) 


The  factor  of  1/2  is  (2.6)  and  in  (2.8)  avoids  double  counting  of  the 
cooperative  terms.  The  assumption  implicit  in  (2.8)  is  that,  having  singled  out 
the  ith  beam  particle,  the  remainder  of  the  beam  is  treated  as  a  contimum  of 
charge,  i.e.  a  charged  fluid  of  density  p(z,b). 


In  order  to  proceed  further,  the  form  of  p(z,b)  is  specified  as 


p(z,b) 


n  e 
B 


e(z+i.1)e(L2-z) 


(2.9) 
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Thus  a  beam  having  a  step  function  [ © C z ) ]  head  »od  tail  and  a  Gaussian  radial 

profile  will  be  considered  in  order  to  illustrate  the  properties  of  the 

beam -aensity  effect.  Lj  is  the  distance  from  the  particle  of  interest  to  the 

beam's  head  and  is  the  distance  to  the  tail.  The  particle  is  assumed  to  be 

on  axis.  Setting  Z.*Z.  and  using  C 2 . 1 )  and  (2.9)  in  (2.8),  W  becomes 
11  B 


Wg  *  ^e2 

If 


x  (1-B2e(v))  K0 


(  bw 

W' 


(l-fi2e  (v)) 


(2.10) 


Next,  the  order  of  the  frequency  and  space  integrals  is  interchanged.  This  is 
permissable  since  the  integrand  is  uniformly  convergent  for  finite  z  and  h.  The 
z-integration  yields 


(2.11  ) 


The  b-integration 


.  23 
is 


-(b/bQ)2 

db  b  e 


K0<b/B) 


b2  (b0/2B)2 
_2  e  Ei 

4  1 


(2.12) 


where  Ej(Z)  is  the  exponential  integral  function,  and 


.  w  v  1/2 

B-1  =  _£_  ( l-(32c )  .  (2.13) 
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The  correction  to  the  integral  for  b  <  a.  Equation  (2.1),  haa  been  neglected 
aince  it  is  small.  The  8-integration  gives  a  factor  of  2i . 


Using  C2.ll)— (2.13)  in  (2.10),  Wg  becomes 


W 

B 


n-b^  2  a)  f 

JL°  e  — £  Re  J 
v 


dv_l__(  1  -B2e  (v  )  ) 
c  (v) 


/b  u  v \ 


(1-S2c ) 


[{YYf 


X 


_e 


iu  L  v 


v 


iu  L  v 
_£_! 
v 


-iw  L  v 
.  P..  1 


(2.14) 


THE  FREQUENCY  INTEGRATION 

At  this  point  the  form  of  the  dielectric  function  must  be  specified.  The 
Lorentzian  dielectric  function  of  the  form 

N 

e(v)  -  1  +  Y  _ i _  (2.15) 

v4— 2in  v-v4- 
i  i 

is  chosen.  In  (2.15),  the  are  the  eigenfrequencies  of  harmonic  oscillators 
which  are  given  values  to  correspond  to  atomic  transitions,  and  the  n£  are  the 
damping  coefficients,  or  the  half  widths  of  the  atomic  transitions.  The 
oscillator  strengths  f^  are  chosen  to  equal  the  fraction  of  an  atom's  electrons 
in  a  subshell.  For  conduction  electrons  *  0  and  2n  .  is  the  plasmon  line 
width  in  solids  or  the  collision  frequency  in  a  plasma,  all  in  units  of  In 

short,  the  parameters  in  c(v)  are  chosen  phenomenologically  to  represent  the 
dielectric  properties  of  various  media. 
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Following  the  method  developed  by  Sternheimer  in  hi#  calculations  of 
the  Fermi  density  effect, 


e-1  (v)  fs  1  - 


£*1  i?-2in.v-v2 

l  l 


(2.16) 


2  2  2 

where  l.  "  p  vt  +  f ..  The  Sternheimer  factor  p  is  chosen  so  that  the 

value  of  the  Bethe  logarithm,  lnl,  obtained  in  non-relativistic  experiments,  is 

reproduced  by 


w 

filndi' 


[Kujp)  3  lnl. 


(2.17) 


Defining 


,  2  21/2  .  .  2  21/2 

*j  ■  ~inj +  (1rnj)  ♦  Bj  *  -inj  -  (1j  • n * 1  j>n j * 


,  2  21/2  ,  ,  2  21/2, 

*  -i[nj-(nj-ij)  ],  8j  *  -i[nj  ♦  (nj-tj)  ],  ijinj, 


Equation  (2.16)  can  be  rewritten  as 


.-!(»>  ■  1  *  T, 

J-l  J  J 


therefore  c  (v)  has  poles  at  v-a^  and  v*Sj  in  the  lower  half  complex 
v-plane. 


(2.18a) 


(2.18b) 


(2.19) 
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In  order  Co  accomplish  Che  frequency  integral  of  EqueCion  (2.14),  two 
contour  integrals  will  be  performed,  one  in  Che  firaC  quadrant  and  one  in  the 
fourth  quadrant  of  the  complex  v-plane.  Therefore,  the  following  integrals 
are  defined: 


I(C+)  2  y  dv  I  (l-fi2e ) 


Av2(1-62 


c )  r  iBiv  iB2vi 

EjIA  v2(l-fi2e)j  je  ♦  e  J 


(2.20) 


C+ 


Av2(l-fi2c)  .  .  r~lBlV  “lB2vl 

(C-)  2  -  f  dv  I  (1-S2c)  e  E^A  v2(1-B2c)J  |e  ♦  e  J 


C_ 


where 


*  ('1^)  • 


(2.22) 


w  L.  i»  L_ 

Bj  5  -P-1  and  B2  5  _E_£  . 

v  v 


(2.23) 


Refering  to  Figure  1,  the  integral  around  the  first  quadrant  is  I(C+)*0  so  chat 


i(c1)  »  -r(c3)  -  r(c2). 


(2.24) 


N 

The  integral  around  Che  fourth  quadrant  is  I(C_)  *  -2ti  /  „  Ri,  where  the 

j-1 

Rj  are  the  residues  of  the  poles  of  t”^(v).  Therefore, 


I(C8) 


-  I(C4)  -  I(C5)  -  I(C6)  -  I(C7). 


(2.25) 


9 


Using  (2.24)  «nd  (2.25)  in  (2.14),  one  nay  write  Wg  aa 


»B 


n-b2  «  w 

JL £  e2  _E  Re  [I(Ci) 
2  v 


♦  I(C8)]. 


(2.26) 


The  contour  integrals  required  to  evaluate  I(C^)  and  I(Cg)  will  be 
discussed  next.  As  the  radii  of  the  quarter  circles  for  and  become 
infinite,  KCj)-^  and  I(C^)-*0.  Along  Cg,  v*iy  and  all  the  factors  of  the 
integrand  are  real  except  dv*idy  and  E^(Z),  where 


Z  s  Av2(l-fi2c). 


(2.27) 


2  24 

There  are  zeros  of  1-B  c  on  the  imaginary  axis  ,  so  that 


1  -  B  c(+il+)  5  0 


(2.28) 


defines  the  positions  v«+il+  and  v**-ii_  of  the  zeroes.  Furthermore,  one  has 


and 


1-B  c(+iy)>0  y>l+ 


l-B  c(+iy)<0  y<l+. 


(2.29a) 


(2.29b) 


which  gives  Z<0  for  y>f+  and  Z>0  for  y<l  +  •  Hence, 


{“Ei(lZl)Tiw  ,  y>i+ 
EjdZl)  y<l  + 


(2.30) 


Taking  (2.30)  into  account, 


(2.31) 


I(C3>  ...  f  dy-^  (,-»*««,>)  r‘y2<1-e2*  ’(.'V  .  .’V, ; 


and  similarly 


m 

/,  /  ~  \  -Ay2(l-B2e)  -B.y 

dy  /l-B^c(-iy)j  e  [e  *  +  e 


-V, 


(2.32) 


It  can  be  shown  that  0<  !+<_  l_<  (l-E2)”^2. 

The  expression  l-82e  can  be  rewritten  as 

r  2N'2  /  N  ,  i 

1-S2e  ■  (1-B2)  l(w-il  +  )(v+il-)  IT  (v-Zq)/  IT  (v+Uj)(v-Uj)l  , 

L  q-1  /  j-1  J 

where  the  factors  in  the  numerator  give  the  positions  of  the  zeroes  of 
2 

1-8  c,  while  the  poles  of  c  are  given  by  the  U/s.  With  the  exception  of 
it+,  all  the  zeroes  occur  in  the  lower  half  plane,  as  do  the  poles  which  are  at 


(2.33) 


(-ini  +  (v2-n2)^ 
v  -  0:  -  ;  j  j 

(-i[n,  -  (n2-v2)] 


(-inj  +  (v2-^2)1^2  vj>nj 
^2]  v?<th 


(2.34) 


and 


!-V*  -  - 

-ifn.  ♦ 


J  J 


-V*  »  -in  -  (v2  -  n2)^2  v  ;>n  i 

4  j  j  j  J  J 

(n?  -  v2)1/2] 
j  J  J 


(2.33) 


The  complex  values  of  Z.  are  obtained  numerically  when  N  is  greater  than  one. 
For  each  pole  U.  in  the  fourth  quadrant,  there  will  be  an  associated  zero  Z. 
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•uch  that  lia  Zi  ■  0:. 

8*0 

In  Figure  1,  the  contour  Cg  encloses  the  pair  of  points  Uj  and  Zj.  Since 
E^(Z)  is  defined  with  a  branch  cut  from  aero  to  a  cut  between  Z^  and  0^  has 
been  introduced  in  Figure  1  such  that  Cg  does  not  cross  the  cut.  The  function 
E^(Z)  on  either  side  of  the  cut  differs  by  2»i,  where  Z  is  defined  by  (2.27). 
Using  this  result  in  (2.21),  the  contour  integral  I(C.)  can  be  expressed  as 


1(C6)  -  2wi 


®j 

i  J  dv  i  (1-B2c )e 


Av2(l-a2c),  -iB.v  -iB  v 
[e  1  +  e  2  ] 


(2.34) 


along  a  path 


w  -  Zi  +  reia. 


where 


6  “  tan' 


.j  r-i«(urz.)  ] 

[  Re  COj-ZjjJ 


(2.35) 


(2.36) 


For  each  pair  (Zj,U.)  in  the  fourth  quadrant  there  will  be  a  contribution  like 

I(c6>. 


If  the  dielectric  function  contains  conduction  electrons,  there  is  a  pole 
of  c  on  the  imaginary  axis  at  v*-2inc,  where  2nc«a^  is  the  collision  frequency 
or  the  plasnon  width,  as  appropriate.  In  this  case  Z  <  0  for  y  <_  2nc  (see 
Equations  (2.27)  and  (2.30))  and  again 
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Sj(Z)  •  -EiOzO  +  is,  y  £  2 nc; 


(2.37) 


therefore. 


I(C7) 


/  ,  ,  -Ay2(l-B2c)  -By  -B,y 

/  dy  (l-B2e(-iy))c  fe  1  +  e  2  ] 

J 


(2.38) 


Note,  for  the  case  in  which  c  contains  conduction  electrons  lim  t.  “  2nc< 

6*0 


The  contribution  to  (2.25)  from  the  residues  Rj  of  the  poles  of  c  at  v  ■ 
o,  (see  (2.19))  can  be  written  by  inspection  of  Equation  (2.21):  if 


-f.  A a?  2  -iB,a .  -iB,a  . 

R;  ■  _ 1L-  e  J  Ei(Aa:)[e  1  J  +  e  2  3] 

°rB  j 


(2.39) 


if 


-f.  -Al  B .  1 2  - 

Us  *  _ — _J e  J  Ei(AiBii2)[e  ♦  e  ] 

J  2(aj-Bj7  J 


-Ala .  1 2  ,  “Bjlajl  "B2laj' 

-e  J  Ei(Alajl2)[e  ♦  e  ] 


(2.40) 


In  the  latter  case  the  poles  lie  on  the  imaginary  axis  (Equation  2.18b). 

All  of  the  quantities  required  to  evaluate  I(C^)  and  I(Cg)  appearing  in 

Wg,  Equation  (2.26),  have  now  been  specified.  I(Cg),  KCj),  I(C^)  and  I(C^) 

have  been  expressed  as  integrals  which  must  be  done  numerically.  In  the 

nonrelativistic  regime,  an  analytic  approximation  to  Wg  has  been  given 
21 

previously,  which  agrees  well  with  the  corresponding  numerically  integrated 


result. 
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CHAPTER  3 

NUMERICAL  RESULTS  FOR  A  RELATIVISTIC  BEAM  IN  NITROGEN 

The  expression  for  W^,  Equation  (2.26)  has  been  evaluated  numerically  for 

the  case  of  an  electron  beam  pulse  with  an  energy  of  50.59  MeV  (fl*0.99995,  y  * 

100.0).  The  Gaussian  radial  profile  parameter  bQ  of  Equation  (2.9)  was  taken 

13  -3 

to  be  0*1  cm  and  the  beam  density  n&  was  6.623  x  10  cm  ,  which  corresponds  to 

5  2 

a  total  current  of  10  KA,  or  a  current  density  of  3.18  x  10  A/ cm  .  The  pulse 
length  was  1  meter. 


The  medium  is  characterized  by  the  dielectric  funccion  e(\>)  with  the 

parameters  discussed  in  connection  with  Equations  (2.15)-(2.17).  Table  1 

contains  the  values  used  in  these  expressions  for  e(v)  and  c  ^(v).  The  first 

three  eigenfrequencies  were  taken  from  Reference  25.  A  more  recent  set  has  been 
•  26 

given  by  Steraheimer  and  Peierls,  but  the  ones  of  Reference  25  are  adequate 

for  the  present  exploratory  purpose.  The  R-shell  line  width  was  obtained  from 
.  27 

the  compilation  of  widths  made  by  Krause,  while  the  L^  and  L^  widths 
were  estimated  from  the  formula  for  collisional  line  widths  of  Weisskopf  as 
applied  in  Reference  28: 


"finj 


f.  (fi  «  )2 

-1 _ E _ 

8  15  Wj 


teV] 


or 


f . 

nj  *  — ,  units  of  Up. 


(3.1) 


(3.2) 
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TABLE  1.  INPOT  PARAMETERS  TO  CHARACTERIZE  THE  N2  MEDIUM 

See  Equation  (2.15)  for  definitions.  The  frequencies  Vj  and  Hj  are  in 

units  of  u>,  where  ‘Rw  *  0.7208  ev.  The  gas  density  is  1.1652  x  10~3 
P  P 

g/cm  .  The  parameter  p  *  1.23,  see  Equations  (2.16)  and  (2.17). 


j 

f. 

v  . 

n  . 

J 

3 

J 

1 

0.2587 

572.4 

4.9  x  10“2 

2 

0.2587 

66.1 

4.4  x  10-4 

3 

0.4286 

-5 

43.4 

1.0  x  10-3 

4 

2.659x10 

0.0 

3.63  x  10~ 
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The  value  for  f^  was  obtained  by  assuming  that  the  M2  medium  was  weakly 

16—3 

ionized  and  had  a  free  electron  density  of  ng  *  10  cm  ,  corresponding  to 
XT  £$0.7  eV.  The  parameter  in  units  of  is  half  the  collision 

frequency  which  can  be  obtained  from  the  momentum  transfer  cross 

29  -15 

section,  a  ■  1,1  x  10 
m 


v 


c 


■  1.0  x  10*2  sec”* 


) 


where  n  is  the  gas  density  (see  Table  I). 


(3.3) 


According  to  the  discussion  below  Equation  (2.16), 


1/2  -3 

l,  -  V,  -  5.156  x  10  , 

4  4 


which  corresponds  to  a  plasma  frequency  of 


*Rw.  -  t.-fiw  «  3.716  x  10-3  eV 
4  4  p 


for  the  free  electron  component.  Note  that  the  ratio 


(3.4) 


(3.5) 


n4/*4  *  0.7, 


(3.6) 


so  that  the  life  time  is  of  the  order  of  the  period  of  oscillation  for  this 
choice  of  parameters. 


17 


9 


t 
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Figure  2  illustrates  the  results  of  the  calculation  of  Wfil  the 
beam-density  effect  on  ionization  energy  loss  for  the  particular  set  of 
parameters  discussed  above.  The  ratio  of  the  beam  to  the  single  particle  term, 
Wg/S,  [see  Equations  (2.3)  and  (2.7)}  is  plotted  versus  the  distance  from  the 
beam-front  in  centimeters.  The  curves  shown  are  for  a  particle  on  the  beam 

axis.  The  single-particle  energy  loss  S  (Equation  (2.3))  was  calculated  using 

24  -2 

Sternheimer *s  method,  with  the  result  that  S“2.48  MeV/g-cm  .  The  curves  are 

symmetric  with  respect  to  the  midpoint  of  the  pulse  at  50  cm.  The  solid  curve 

16  —3 

was  generated  using  the  parameter  values  given  in  Table  1  with  ng  =  10  cm  , 

17  -3 

while  the  dashed  curve  results  when  n  ■  10  cm  and  the  dash-dot  curve  was 

e 

15  -3 

obtained  with  n  *  10  cm  .  As  above,  n  is  the  number  of  free  electrons 
e  e 

assumed  to  be  present  in  the  medium  through  which  the  pulse  passes  as  it 
undergoes  energy  loss. 

All  three  curves  exhibit  a  large  beam-density  effect  within  the  first 
millimeter  from  the  beam  front.  At  the  midpoint  of  the  pulse  the  three  cases 

yield  a  ratio  Wg/S  of  between  one  and  two,  which  is  still  quite  significant. 

.  .  17  -3 

The  oscillatory  behavior  which  occurs  when  nfi  *  10  cm  comes  from  the 

residues  of  the  poles  of  e  *  given  by  (2.39).  The  residue  terms  decay 

exponentially  and  begin  to  be  dominated  by  the  nonresonant  terms  at  about  0.05 

cm  from  the  beam  front.  When  n  *  10^  and  10^  cm  \  the  nonresonant 

e 

contribution  dominates  the  resonant  part  and  the  oscillatory  behavior  disappears. 

For  a  particle  at  the  beam  front,  *  0  and  is  equal  to  the  pulse 
length.  Equation  (2.9).  For  a  semi-infinite  pulse,  then  if  L^O  one 

obtains  Wg"0  (see  Equation  (2.14)).  Even  for  the  one  meter  pulse  studied, 

Wg/S  ~  10  *  at  Lj^O,  Lj*lm.  This  behavior  provides  a  useful  check  on 
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the  numerical  integrations,  since  the  nonresonant  terms  calculated  numerically 

must  cancel  exactly  the  analytical  resonant  terms  at  if  Lj*  ®.  This  test 

.  .  15  -3 

was  particularly  useful  in  checking  the  case  of  ne=10  cm  since  then  ^ 

and  the  poles  of  £  ^(v)  at  v=a^  and  v=8^  (see  Equation  2.18b)  lie  on  the 
imaginary  axis.  Then  a  principal  value  integral  is  required  in  order  to 
evaluate  I(C^),  Equation  (2.38),  and  the  results  are  very  sensitive  to  the 
numerical  procedure  used. 

The  results  obtained  for  the  individual  contours  of  Figure  1  showed  that 

the  contributions  from  1(0^),  I(C^),  and  l(Cg)  were  negligible  in  this  example. 

I(C^),  however,  resulted  in  the  large  beam-density  contributions  illustrated  in 

Figure  2.  When  0<<1,  i+*0  and  i_*2n^  and  there  will  be  a  small  contribution 

from  I(Cj);  also  as  0*0  Zj*Uj  and  l(Cg)*0.  It  should  be  pointed  out  that  the 

parameters  used  to  model  the  dielectric  function  e (v )  given  in  Table  1  are  based 

on  the  ionization  potentials  of  the  subshells  of  nitrogen  and  are  not  suitable 

for  the  investigation  of  Cerenkov  radiation  by  beams;  therefore,  one  does  not 

expect  to  obtain  any  beam-density  effect  on  Cerenkov  radiation  in  the  present 

28 

example.  By  a  suitable  modification  of  e(v),  the  Cerenkov  part  of  energy  loss 
can  easily  be  incorporated  into  this  treatment. 

The  contributions  to  Wg  from  the  contour  integrals  of  Figure  1  all  contain 

both  longitudinal  and  transverse  field  components.  The  transverse  components 

2  2 

come  from  the  singularities  associated  with  E^(Av  (1-B  e)),  Equation  (2.21),  as 

2  ... 

1-0  c*0  and  £*®.  In  the  limit  0*0  only  longitudinal  contributions  remain.  In 


Reference  21  it  was  shown  that  the  longitudinal  contribution  from  the  contour 
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along  the  imaginary  axis  was  of  a  long-range,  dipole-like  form  for  the 
cooperative  loss  between  two  particles.  It  is  this  long  range,  nonresonant 
contribution  which  is  dominant  in  the  present  calculation.  The  resonant  part, 
coming  from  the  poles  of  e  *(v)  Equation  (2.39)  is  purely  longitudinal  and 
is  a  relatively  short  ranged  term. 
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CHAPTER  4 
DISCUSSION 


The  evaluation  of  W^,  Equation  (2.26),  and  the  results  of  the  beam-density 
effect  calculation  presented  in  Chapter  3  contain  no  approximation  other  than 
those  implicit  in  the  models  chosen  to  represent  the  medium  and  the  beam.  The 
medium  was  characterized  by  a  dielectric  function  based  on  the  assumption  of 
linear  response  to  the  local  IT-  field.  The  validity  of  this  assumption  may  be 
questionable  near  the  beam-front  where  the  fields  have  large,  rapid  variations. 
The  derivation  of  Equation  (2.26)  was  based  on  a  spacially  non-dispersive 
dielectric  function,  but  a  minimum  impact  parameter  cut  off  was  used  (see 
Equation  (2.1))  to  separate  the  region  of  distant  collisions  with  small  momentum 
transfers  from  the  region  of  high  momentum  collisions.  This  procedure  has  been 
found  to  yield  good  agreement  with  two-particle  vicinage  function  calculations 
using  a  simple  momentum  dependent  dielectric  function  in  the  nonrelativistic 
regime  (see  References  9-12  and  19).  After  integration  over  the  beam,  the 
effect  of  a  non-zero  minimum  impact  parameter  was  negligible. 

The  medium  was  assumed  to  remain  in  its  initial  state  during  the 
interaction  with  the  pulse.  Therefore,  the  increase  in  medium  conductivity 
produced  by  the  deposited  beam  energy  has  not  been  incorporated  in  a  self- 
consistent  way  with  the  energy-deposition  calculation.  The  variation  of  ng, 

3 

the  number  of  free  electrons  per  cm  ,  was  given  in  Figure  2  to  illustrate 
roughly  the  possible  effect  of  changes  in  conductivity.  This  first  calculation 
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is  intended  primarily  to  explore  the  basic  properties  and  magnitude  of  the 
beam-density  effect  and  its  relation  to  other  phenomena  rather  than  to  be  a 
definitive  self-consistent  treatment. 

The  beam  pulse  was  modeled  by  a  charged  fluid  with  a  step  function  head  and 
tail  and  a  Gaussian  radial  profile.  The  abrupt  rise  of  current  in  this  model  is 
somewhat  unrealistic.  For  a  more  gradual  current  rise  and  fall,  the  values  of 
Wg  near  the  beam  front  and  tail  may  be  considerably  less  than  those  found  in 
the  present  calculation.  The  question  of  how  well  a  fluid  model  simulates  the 
more  realistic  case  of  a  system  of  discrete  beam  particles  which  are  distributed 
according  to  an  appropriate  correlation  function  requires  further  consideration 
also.  Note  that  the  starting  point  of  the  present  treatment  was  two  discrete 
beam  particles  separated  by  an  arbitrary  position  vector. 

Setting  aside  the  caveats  discussed  above.  Figure  2  shows  that  the 
particles  in  a  1  meter  pulse  experience  energy  loss  at  a  rate  which  is  greater 
than  twice  that  of  a  single-particle.  Near  the  head  and  tail,  the  loss  rate  for 
a  50  MeV,  10KA  beam  of  0.1  cm  radius  can  range  over  one  or  two  orders  of 
magnitude,  depending  on  the  ionisation  state  of  nitrogen  at  1  atmosphere. 

Longer  pulses  will  have  a  smaller  proportion  of  their  electrons  which  experience 
a  large  beam-density  loss.  The  magnitude  of  this  loss  per  particle  is  directly 
proportional  to  the  current  density.  Further  studies  of  the  dependence  of  Wg 
on  medium  and  beam  parameters  are  underway. 

In  order  to  relate  the  beam-density  effect  as  presented  in  the  present 

paper  to  the  single-particle  ionization  energy  loss,  and  at  the  same  time  to  the 

1 A  X3 

phenomena  of  ohmic  loss  and  return  currents1  ’  in  plasmas,  a  brief  discussion 
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of  ch«  various  approaches  to  energy  loss  will  be  useful.  The  ionization  energy 
loss  of  a  single  particle  may  be  considered  from  any  of  three  equivalent  points 
of  view:  a)  the  work  done  on  the  medium  by  the  particle's  7-field,  the  Bohr 
viewpoint,  b)  the  work  done  on  the  particle  by  the  field  induced  in  the 
medium,  the  Landau  viewpoint,  and  c)  the  Poynting  flux  through  the  surface  of  a 
cylinder  surrounding  the  particle's  trajectory,  the  Fermi  viewpoint.  The 
equivalence  of  these  three  approaches  is  well  known  and  can  be  proven  using  the 
Poynting  theorem.  In  the  case  of  beams  of  particles,  a  similar  equivalence  has 
been  derived  in  Reference  19.  There  it  was  shown  that  the  rate  of  change  of 
mechanical  energy  of  a  beam  caused  by  the  induced  field  of  the  medium  is 


-  /7  •  (7-7)d3x  -  -4w  }  Zievi  *  T, 

dt  J  u 


(4.1) 


where  P  is  the  total  induced  polarization  field  caused  by  all  beam  particles; 
4*P**1?-If;  and  the  beam  current  is  7*  ^  *  Z ^ev^S  Cr-v^t-R^)  ♦  This  expression 
excludes  internal  transformation  of  beam  field  energy  to  beam  mechanical  energy, 
which  would  occur  even  in  the  absence  of  a  medium.  In  terms  of  the  induced 


polarization  field  associated  with  each  beam  particle,  P®  /  ] 
stay  write 


and  one 


5T  "  _4w  jL  Zi*Vi  *  * 


(4.2) 


In  this  expression  the  terms  with  iaj  are  the  usual  single-particle  terms  and 
those  with  i*j  are  the  cooperative  beam  terms.  For  the  ith  particle,  one  has 
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-  -4*ZieVi  •  7i  -4*ZieVi  •  )  i  Pj  ,  j  *»i.  (4.3) 

J=1 

If  all  Che  bean  particles  are  assumed  Co  have  equal  velocities,  v,  Chen 
V"ldEHi  /dt  corresponds  Co  EquaCion  (2.7). 

Focusing  aCCention  on  Che  beam,  one  sees  Chat  (4.1)  or  (4.3)  has  Che 

general  form  of  ohmic  loss  caused  by  Che  inCeraction  of  Che  induced  medium  field 

with  Che  beam.  This  is  essentially  Che  Landau  viewpoint  discussed  above.  If, 

instead,  Che  response  of  Che  medium  Co  Che  beam's  fields  is  considered,  Chen  one 

has  Che  Bohr  viewpoint.  Since  Che  dieleccric  function  used  in  Che  present 

calculation  includes  Che  contribution  from  Che  conduction  electrons,  the  Bohr 

viewpoint  will  encompass  the  phenomenon  of  the  return  current  induced  by  the 

14 

beam,  as  well  as  Che  polarization  currents  of  the  bound  electrons.  These 
currents  are  just  dP*/dt. 

In  Che  presenC  calculation  it  was  found  Chat  the  dominant  beam-density 
effecC  comes  from  the  nonresonant,  low  frequency  contribution  to  the  frequency 
integral  of  W .  This  part  of  the  effect  may  be  identified  with  ohmic  loss  in 

D 

a  plasma.  However,  both  the  resonant  and  nonresonant  terms  are  required  in 

order  to  obtain  the  proper  behavior  of  Wg  at  the  beam  front  and  at  a  point 

infinitely  far  from  the  head  or  tail  of  an  infinite  beam.  The  relation  of 

V(z,b),  Equation  (2.1),  to  the  nonrelativistic  two-particle  vicinage 
7  19 

functions  ’  has  already  been  discussed  in  Chapter  2.  Formally,  W(z,b)  also 

.  .  20  v 
contains  the  Fermi  density  effect,  including  the  Cerenkov  loss  term.  The 

integral  of  W(z,b)  over  a  beam  further  incorporates  the  ohmic  loss  term  along 

with  these  one-  and  two-  particle  energy  loss  phenomena.  Therefore,  from  a 

formal  point  of  view  Che  present  treatment  includes  all  the  various  aspects  of 

collisional  energy  loss  by  charged  particle  beams  in  a  unified  way. 


HSWC  TR  83-348 


REFERENCES 

1.  Toepfer,  Alan  J.,  "Particle  Beam  Fusion, "  Advances  in  Electronics  and 
Electron  Physics.  Vol.  53,  1980,  p.l. 

2.  Imasaki,  R.,  et  al.,  "Evidence  of  Anomalous  Interaction  Between  a 
Relativistic  Electron  Beam  and  Solid  Target,"  Physical  Review  Letters,  Vol. 
43,  1979,  p.  1937. 

3.  Imasaki,  K.,  et  al.,  "Energy  Absorption  and  Transport  in  Layered  Targets 
Irradiated  by  a  Relativistic  Electron  Beam,"  Applied  Physics  Letters.  Vol. 
37,  No.  6,  1980,  p.  533. 

4.  Young,  F.  C.  et  al.,  "Measurements  of  Enhanced  Stopping  of  1-MeV  Deuterons 
in  Target-Ablation  Plasmas,"  Physical  Review  Letters,  Vol.  49,  1982,  p.  549. 

5.  Kiselev,  V.  A.,  et  al.,  "Interaction  of  a  Relativistic  Monochromatic 
Electron  Beam  with  Plasma  Produced  in  the  Atmosphere,"  Soviet  Physics  JETP 
Letters,  Vol.  29,  No.  12,  1979,  p.  700. 

6.  McCorkle,  R.  A.,  and  lafrate,  G.  J.,  "Beam-Density  Effect  on  the  Stopping 
of  Fast  Charged  Particles  in  Matter,"  Physical  Review  Letters.  Vol.  39,  No. 
20,  1977,  p.  1263. 


I'l  I  iv 


27 


NSWC  TS  83-348 


REFERENCES  (Cont.) 

7*  Brandt,  Werner,  Ratkowski,  Anthony,  and  Ritchie,  R.  H. ,  "Energy  Loss  of 

Swift  Proton  Clusters  in  SoLids,"  Physical  Review  Letters,  Vol.  33,  No.  22, 
1974,  p.  1325. 

8.  Brandt,  Werner  and  Ritchie,  R.  H. ,  "Penetration  of  Swift  Ion  Clusters 
Through  Solids,"  Nuclear  Instruments  and  Methods,  Vol.  132,  1976,  p.  43. 

9.  Basbas,  George,  and  Ritchie,  R.  H. ,  "Vicinage  Effects  in  Ion-Cluster 
Collisions  with  Condensed  Matter  and  with  Single  Atoms,"  Physical  Review  A, 
Vol.  25,  No.  4,  1982,  p.  1943. 

10.  Neelavathi,  V.  N. ,  Ritchie,  R.  H. ,  and  Brandt,  Werner,  "Bound  Electron 
States  in  the  Wake  of  Swift  Ions  in  Solids,"  Physical  Review  Letters,  Vol. 
33,  No.  5,  1974,  p.  302. 

11.  Ritchie,  R.  H. ,  Brandt,  Werner,  and  Echenique,  P.  M. ,  "Wake  Potential  of 
Swift  Ions  in  Solids,"  Physical  Review  B,  Vol.  14,  No.  11,  1976,  p.  4808. 

12.  Mazarro,  A.,  Echenique,  P.  M. ,  and  Ritchie,  R.  H. ,  "Charge-Particle  Wake  in 
the  Random-Phase  Approximation,"  Physical  Review  B,  Vol.  27,  No.  7,  1983, 
p.  4117. 

13.  McCorkle,  R.  A.,  and  Cox,  J.  L.  Jr.,  "Beam-Plasma  Electrodynamics,"  Physics 
Letters,  Vol.  83A,  No.  9. 


28 


NS WC  TR  83*348 


REFERENCES  (Cont.) 

14.  Roberta,  T.  G.,  end  Bennet,  W.  H.,  "The  Pinch  Effect  in  Pulsed  Streams  at 
Relativistic  Energies,"  Plasma  Physics.  Vol.  10,  1968,  p.  381. 

* 

15.  Cox,  James  L. ,  Jr.  and  Bennett,  Hillard  H.,  "Reverse  Current  Induced  by 
Injection  of  a  Relativistic  Electron  Beam  into  a  Pinched  Plasma,"  Physics 
of  Pluids.  Vol.  13,  No.  1,  1970,  p.  182. 

16.  Hasmier,  D.  A.,  and  Rostocker,  N. ,  "Propagation  of  High  Current  Relativistic 
Electron  Beams,"  Physics  of  Fluids,  Vol  13,  No.  7,  1970,  p.  1831. 

17.  Lovelace,  R.  V.,  and  Sudan,  R.  N. ,  "Plasma  Heating  by  High-Current 
Relativistic  Electron  Beams,"  Physical  Review  Letters,  Vol.  27,  No.  19, 
1971,  p.  1256. 

18.  Wallis,  G.  et  al.,  "Injection  of  High-Current  Relativistic  Electron  Beams 
into  Plasma,"  Soviet  Physics-Pspekhi ,  Vol.  17,  No.  4,  1975,  p.  492. 

19.  Rule,  D.  H.,  and  Cha,  M.  H.,  "Collective  Effects  in  Energy  Loss  by 
Relativistic  Clusters  of  Charged  Particles,"  Physical  Review  A,  Vol.  24, 

No.  1,  1981,  p.  55. 

20.  Rule,  D.  W.,  "Collective  Effect  in  Energy  Deposition  by  Beams  of 

•  Relativistic  Charged  Particles,"  in  High-Power  Beams  81,  Proceedings  of  the 

Pourth  International  Topical  Conference  on  High-Power  Electron  and  Ion-Beam 
Research  and  Technology-Palaiseau,  29  June  1981,  pp.  481-488. 

29 

L _ 


4 


NSWC  TR  83-348 


REFERENCES  (Cont.) 

21.  Rule,  D.  W.,  and  Crawford,  0.  H. ,  The  Nature  of  the  Beam-Density  Effect  on 
Energy  Loss  bv  Nonrelativistic  Charged  Particle  Beams.  NSWC  TR  83-346, 

June  1983. 

22.  Jackson,  J.  D.,  Classical  Electrodynamics,  2nd  ed.  (John  Wiley  and  Sons, 
Inc.,  New  York,  1975)  pp.  618-653. 

23.  Gradshteyn,  I.  S.,  and  Ryzhik,  I.  M. ,  Tables  of  Integrals,  Series,  and 
Products  (Academic  Press,  New  York,  1980),  p.  717. 

24.  Sternheimer,  R.  M.,  "The  Density  Effect  for  the  Ionization  Loss  in  Various 
Materials,"  Physical  Review,  Vol.  88,  No.  4,  1952,  p.  851. 

25.  Sternheimer,  R.  M.,  "Density  Effect  for  the  Ionization  Loss  in  Various 
Materials,"  Physical  Review,  Vol.  103,  No.  3,  1956,  p.  511. 

26.  Sternheimer,  R.  M. ,  and  Peierls,  R.  F.,"  General  Expression  for  the  Density 
Effect  for  the  Ionization  Loss  of  Charged  Particles,"  Physical  Review  B, 
Vol.  3,  No.  11,  1971,  p.  3681. 

27.  Krause,  M.  0.,  "Atomic  Radiative  and  Radiationless  Yields  for  K  and  L 
Shells,"  Journal  of  Physical  and  Chemical  Reference  Data,  Vol.  8,  1979, 
p.  307. 

.  Sternheimer,  R.  M. ,  "The  Energy  Loss  of  a  Fast  Charged  Particle  by  Cerenkov 

Radiation,"  Physical  Review,  Vol.  91,  No.  2,  1953,  p.  256. 

30 


28 


NSWC  TR  83-348 


REFERENCES  (Cont.) 

Itikava,  Yukikazu,  "Momentum  Transfer  Cross  Sections  for  Electron 
Collisions  with  Atoms  and  Molecules,"  Atomic  Data  and  Nuclear  Data  Tables, 
Vol.  14,  No.  1,  1974,  p.  1. 

Landau,  L.  M.,  and  Lifshitz,  E.  M. ,  Electrodynamics  of  Continuous  Media 
(Addison-Wesley  Publishing  Co.,  Inc.,  Reading,  MA,  i960)  pp.  344-359. 


31/32 


NSWC  TR  83-348 


DISTRIBUTION 


Copi es 

4 

Commanding  Officer  Naval  Research  Lab 
Attn:  Code  4700 

,  (Or.  W.A.  All)  1 

Code  6650 

(Dr.  J.B.  Aviles,  Jr.)  1 
Code  4770 

(Dr.  Richard  Fernsler)  1 


(Dr.  S.  Goldstein)  1 

Code  4763 

(Dr.  Robert  Grelg)  1 

Code  4790 

(Dr.  M.  Lampe)  1 

Code  4740 

(Dr.  C.  W.  Roberson)  1 

Code  6654 

(Dr.  Alex  Stolovy)  1 

(Dr.  F.  C.  Young)  1 

Washington,  D.  C.  20375 

Commander,  Naval  Sea  Systems  Command 
Attn:  (D.L.  Merritt)  1 

11N08  NCI 

Washington,  D.  C.  20362 

Office  of  Naval  Research 
Attn:  (Dr.  B.  R.  Junker)  1 

800  N.  Quincy  St. 

Arlington,  VA  22217 

Commander,  Harry  Diamond  Lab 
Attn:  (Dr.  S.  Grayblll)  1 

2800  Powder  Mill  Road 
Adel  phi,  N)  20783 

U.S.  Army  Ballistic  Research  Lab. 
Attn:  (Dr.  D.  Eccleshall)  1 

Aberdeen  Proving  Ground,  MD  21005 

4FWAL-MLP0 

Attn:  (Dr.  Peter  L.  Land)  1 

Bldg.  451 


Wrlght-Patterson  AFB,  OH  45433 

* 

Director,  Nat‘1.  Bureau  of  Standards 
Center  for  Radiation  Research 
*■  Attn:  (Dr.  M.d.  Berger)  1 

(Dr.  S.  M.  Seltzer)  1 

Washington,  D.  C.  20234 


Copies 

Director,  National  Bureau  of  Standards 
Attn:  (Dr.  Mark  Wilson)  1 

Gaithersburg,  MD  20760 

Physical  Review  Letters 

Attn:  (Dr.  George  Basbas)  1 

Box  1000 

Ridge,  NY  11961 

Lawrence  Livermore  Natl.  Lab. 

Attn:  (Dr.  R.  Briggs)  1 

P.0.  Box  808 
Livermore,  CA  94550 

Lawrence  Livermore  Natl.  Lab. 

Attn:  (Dr.  Thomas  Fessenden)  1 

(Dr.  Simon  S.  Yu)  1 

L321 

University  of  California 
Livermore,  CA  94550 

La  jolla  Institute 

Attn:  (Dr.  K.  Brueckner)  1 

P.0.  Box  1434 

La  Jolla,  CA  92038 

Oak  Ridge  National  Lab. 

Attn:  Chemistry  Division 

(Dr.  Oakley  Crawford)  1 
Health  &  Safety  Res.  Div. 

(Dr.  Rufus  H.  Ritchie)  1 
Oak  Ridge,  TN  37830 

Los  Alamos  National  Laboratory 
Attn:  (Dr.  Harold  Dogllani)  1 

Mall  Sta.  H818,  P.0.  Box  1663 

Los  Alamos,  NM  87545 

Physical  Dynamics,  Inc. 

Attn:  (Dr.  George  Gillespie)  1 

P.  0.  Box  1883 

La  Jolla,  CA  92038 

Laboratory  of  Plasma  Studies 
Attn:  (Prof.  D.  Hammer)  1 

Cornell  University 
Ithaca,  NY  14850 


0) 


9 


NSWC  TR  83-348 


DISTRIBUTION  (Cont.) 


Copies 

Argonne  National  Laboratory 
Attn:  (Dr.  Mitio  Inokuti)  1 

Building  203 
9700  South  Cass  Avenue 
Argonne,  ILL  60439 

Science  Applications,  Inc. 

Attn:  (Dr.  Robert  Johnston)  1 

5  Palo  Alto  Square 
Palo  Alto,  CA  94304 

General  Dynamics  Corp. 

Materials  Technology  Div. 

Attn:  (Dr.  Frank  L.  Madarasz)  1 

Mail  Zone  44-68,  P.  0.  Box  2507 
Pomona,  CA  91769 

Sandia  National  Laboratories 
Attn:  (Or.  E.J.  Mcguire)  1 

Target  Int.  Theory  Div.-  4247 
(Or.  T.  A.  Mehlhorn)  1 

Division  5246 
(Dr.  Bruce  R.  Miller)  1 

Albuquerque,  NM  87115 

Thomas  J.  Watson  Research  Center 
I8M 

Attn:  (Dr.  Richard  McCorkle)  1 

Yorktown  Heights,  N.  Y.  10598 

University  of  California 
Dept,  of  Physics 

Attn:  (Prof.  N.  Rostoker)  1 

Irvine,  CA  92664 


Copies 

Dartmouth  College 

Dept,  of  Physics 


Attn:  (or.  John  E.  Walsh)  1 

Hanover,  NH  03755 

Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314  12 

Internal  Distribution 
R  1 

R04  1 

R40  1 

R401  1 

R13  (J.  Forbes)  1 

R41  (P.  W.  Hesse)  1 

R41  (R.  Cawley)  1 

R41  (M.  H.  Cha)  1 

R41  (H.  C.  Chen)  1 

R41  (  R.  Fiorito)  1 

R41  (0.  F.  Goktepe)  1 

R41  (H.  S.  Uhm)  1 

R34  (DR.  J.  Sharma)  1 

F  1 

F34  (E.  Nolting)  1 

F04  (M.  F.  Rose)  1 

E431  9 

E432  3 

E35  1 


Brookhaven  National  Laboratory 
Dept,  of  Physics 

Attn:  (Dr,  R.M.  Sternheimer  1 

Upton,  NY  11973 

Cornell  University 

LaD.  of  Plasma  Studies 
Attn:  (Prof.  R.  Sudan)  1 

Ithaca,  NY  14850 


(2) 


